Competing ferromagnetic and nematic alignment in self-propelled polar particles 
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We study a Vicsek-style model of self-propelled particles where ferromagnetic and nematic align- 
ment compete in both the usual "metric" version and in the "metric-free" case where a particle 
interacts with its Voronoi neighbors. We show that the phase diagram of this out-of-equilibrium XY 
model is similar to that of its equilibrium counterpart: the properties of the fully-nematic model, 
studied before in [8j, are thus robust to the introduction of a modest bias of interactions towards fer- 
romagnetic alignment. The direct transitions between polar and nematic ordered phases are shown 
to be discontinuous in the metric case, and continuous, belonging to the Ising universality class, in 
the metric-free version. 

PACS numbers: 05.65.+b, 87.18.Gh, 64.60.De 



o 



-a 
c 

o 
o 



00 

cn 

O 

(N 



X 



Collective motion is a topic currently enjoying inter- 
est in various communities Within (statistical) 
physics, the seminal work of Vicsek et al. [5[, followed 
by the remarkable calculation of Toner and Tu has 
offered to view the emergence of collective motion in lead- 
erless groups of identical individuals as the spontaneous 
breaking of rotational invariance. The celebrated Vicsek 
model, which consists of self-propelled particles aligning 
ferromagnetically their orientations with that of their 
neighbors in the presence of noise, was originally pre- 
sented — and rightly so — as an out-of-equilibrium XY 
model where spins are forced to move. As is now well- 
known, the Vicsek model is endowed with properties very 
different from those of the XY model: in two dimensions, 
true long-range polar (orientational) order emerges @ 
from a discontinuous phase transition Q, and the long- 
range correlations and anomalous fluctuations predicted 
by Toner and Tu for the ordered collective motion phase, 
although not observed in the region near onset, are in- 
deed present in a large portion of parameter space. 

Other, Vicsek-style, flocking models have been intro- 
duced which serve as key members of different universal- 
ity classes for "dry active matter" i.e. situations in 
which global momentum is not conserved and hydrody- 
namic interactions play no significant role. A prominent 
case is the "self-propelled rods" model, in which the fer- 
romagnetic interaction of the Vicsek model is replaced by 
nematic alignment in line with the typical outcome 
of inelastic collisions between moving elongated objects. 
Switching from ferromagnetic to nematic symmetry of 
interactions in this other out-of-equilibrium XY model, 
changes the symmetry of the ordered phase (which is 
then nematic), in line with the symmetry change of the 
quasi-ordered phase of the corresponding equilibrium ver- 
sion Q . Numerical results Q suggest that like the origi- 
nal (ferromagnetic) Vicsek model, the nematic order ob- 



(a) 




FIG. 1: (color online) (a) Sketch of two self-propelled needle- 
like objects (black arrows) moving at the same speed in some 
overdamped dynamics, which collide with an incoming an- 
gle of exactly I. The middle of each rod is indicated by 
the thin red line. If the impact point is along the first half 
of the hit needle, (some degree of) polar alignment is ex- 
pected (top panel), whereas anti-alignment will typically oc- 
cur for an impact point at the rear (bottom panel). Fric- 
tion forces may very well, though, lead to polar alignment, 
not anti-alignment, even if the impact point is slightly be- 
yond the first half of the needle. In such a case, the nematic 
symmetry of interactions would be weakly broken and biased 
towards ferromagnetic alignment, (b) Schematic phase dia- 
gram of an equilibrium generalized XY model with Hamilto- 
nian H = — Yl(ij) sc °s[#i — 9j] + (1 — s) cos[2(#i — 9j)] (after 
[l|). Interactions are purely nematic (resp. ferromagnetic) for 
s = (resp. 1). D, P, and N, respectively stand for disorder, 
polar order, and nematic order. 



served is truly long-range, but no Toner- Tu-like calcula- 
tion is available to confirm this at some analytical level. 

One might object that the properties of the "nematic 
Vicsek model" are not robust in the sense that the strict 
nematic symmetry of its interactions may be generically 
broken, albeit weakly, by "friction" effects during colli- 
sions between actual rods (Fig. |TJl). On the other hand, 
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in equilibrium, it is known that the quasi-ordered phase 
of the XY model with nematic interactions resists a mod- 
est amount of ferromagnetic alignment (Fig. [lb) [9] . 

In this Rapid Communication, we introduce and study 
an out-of-equilibrium, Vicsek-style, version of such a 
generalized XY model where ferromagnetic and nematic 
alignment compete. We show that its phase diagram is 
similar to that of its equilibrium counterpart, although 
rendered complicated by density-segregated sub-phases. 
Thus, the fully-nematic Vicsek model is robust to the in- 
troduction of a modest bias of interactions towards ferro- 
magnetic alignment, alleviating the concern raised above 
for collisions of actual rods. We also show that this con- 
clusion holds in the case of "topological" neighbors where 
interactions are not limited to some metric zone but occur 
with those objects defining the first shell of Voronoi cells 
around the considered particle. In both the metric and 
this "metric-free" case, we study the direct polar-nematic 
transition present in phase diagrams like that of Fig. [TJd. 
We provide evidence that it seems to be discontinuous in 
the metric model, but continuous in the metric-free case, 
with critical exponents of the Ising universality class, as 
in equilibrium [9J. We finally discuss the inherent diffi- 
culties in deriving hydrodynamic theories in the case of 
mixed ferromagnetic and nematic interactions. 

Our starting point is a Vicsek-style model with com- 
peting ferromagnetic and nematic interactions: N point 
particles move off-lattice at constant speed Vq on a two 
dimensional L x L torus; particle j is defined by its po- 
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(cos 9j , sin 0j) , the sum is taken over all 
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particles k within unit distance from particle j (includ- 
ing j itself), £ is a white noise uniformly distributed in 
[— §, §] , and the complex stochastic function g s (9, 9') is: 
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with prob. s 
otherwise 
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Note that in the second case g s (9,9') is invariant under 
the transformation 9' —> 9' + tt and thus codes nematic 
alignment, while the first case expresses ferromagnetic 
interaction. For s = 0, this model reduces to model 
studied in 0, while s = 1 is fully equivalent to the stan- 
dard Vicsek model. Thus, s is a key parameter governing 
the relative weight of ferromagnetic interactions which 
comes in addition to the two main parameters, the den- 
sity p = N/L 2 and the noise strength r\. In the following, 
we focus on a low density system at p — | with vo = -|, 
as in [8], and study the (s,f?) parameter plane. [13] 

Systematic scans were performed for different sizes. 
While L = 512 allows for a rough determination of the 
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FIG. 2: (color online) Top panels: Phase diagram of the gen- 
eralized Vicsek model with competing ferromagnetic and ne- 
matic interactions (p = |, L = 1024, Vo = |) (the top right 
panel is a close up of the left panel.). Results are presented 
in function of ^/s for clarity. Symbols located where individ- 
ual runs were performed, coding the resulting observed phase: 
red triangles for polar order, black circles for nematic order, 
and cyan diamonds for disorder. Full/empty symbols code 
for homogeneous/density segregated phases. Bottom panels: 
typical snapshots of the 3 density segregated phases present 
(system size L = 1024) Left: P W avos at s = 0.05, -q = 0.1. Mid- 
dle: Nband at s = 0.01, t) = 0.14. Right: N chaos at s = 0.01, 
r\ = 0.166. The thick magenta arrows indicate the direction 
of motion of particles. 



main features of the phase diagram, one needs larger 
sizes to capture its details, mostly because the segregated 
banded states of the fully nematic (s = 0) model arise 
clearly only for large systems Q . Polar and nematic order 
were characterized by means of the two time-dependent 
global scalar order parameters Pit) = | (exp(i6 , *)) J - 1 (po- 
lar) and Q(t) = |(exp(i26'*)) : ,| (nematic), as well as their 
asymptotic time averages P = (P(t)) t and Q — (Q(t)) t - 
To detect the various density-segregated phases, we re- 
lied on snapshots and movies of coarse-grained density 
and order parameter fields. 

The phase diagram is presented in Fig.[5]for L — 1024, 
a size beyond which its various features do not seem 
to change. Close to the s = axis, one observes that 
the various nematic sub-phases described in [1] for the 
purely nematic case are extended to finite but small s 
values: the spatially-homogeneous nematically-ordered 
phase (denoted Nh om ) can be observed up to ^fs rj 0.17 
for r] « 0.1. The segregated phase with a unique, dense, 
nematically-ordered band occupying a finite fraction of 
space (Nband) extends up to y/s w 0.4 for 77 « 0.13. The 
spectacular spatiotcmporal chaos regime in which thin 
unstable dense bands elongate, twist, break, collide and 
form again (denoted N c h a0 s) is limited to a narrow tongue 
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near r\ = 0.16. As s is increased, the width of this tongue 
decreases, making it difficult to locate numerically; our 
results suggest, though, that it extends all the way to the 
meeting point with the polar order region. 

Polar order can be observed at arbitrary small s, pro- 
vided r\ is small enough; it is delimited by a smooth line 
(like in the equilibrium case, see Fig. [2]) starting at the 
origin and ending at r\ ~ 0.4 for s = 1, the transition 
point of the Vicsek model at this density [ill]. The polar 
order region is itself divided in two, as expected from our 
knowledge of the Vicsek model: in a large, tongue-like re- 
gion bordering the onset of polar order (P W aves in Fig. [2]) , 
the ordered phase consists of trains of solitary traveling 
waves, and is thus different from the Toner- Tu homoge- 
neous phase (Phom) present at lower noise strength. 

The phase diagram of our out-of-equilibrium, Vicsek- 
style, generalized XY model thus possesses the same gen- 
eral structure as its equilibrium counterpart: a small ne- 
matic triangular region is present on the small-s side 
above a continuous line delimiting polar order. This 
general structure is complicated by the presence of the 
various density-segregated ordered subphases. At our 
numerical resolution the line dividing the polar region 
(red dashed line in Fig. [5]) and that dividing the nematic 
region (black dashed line) seem to meet the border of 
the polar region (blue solid line) at the same point flij ]. 
There are thus two different ways of transitioning directly 
from polar to nematic order, as opposed to only one at 
equilibrium: at small s values, the P-N transition occurs 
between spatially-homogeneous phases, while at interme- 
diate s values, it links the two density segregated phases 
P wavos and Nband in Fig. El 

We now turn our attention to the nature of the main 
transitions (disorder/nematic, disorder/polar, and po- 
lar/nematic). Let us first recall that the ordered phases 
observed all seem numerically to show true long-range 
order, even the nematic phases Nhom and Nband @- In 
these last cases, however, as argued in Q, one cannot 
exclude the possibility that nematic order is only quasi- 
long-range at very large scales. In the following, we as- 
sume true long-range order, in agreement with the nu- 
merical results. 

The transition from disorder to nematic order actually 
occurs between the chaotic N c h a os phase and the ordered, 
segregated Nband phase. As reported in |8j, it is deter- 
mined by the long-wavelength instability of the band, 
which leads to the N c haos phase. This phase being disor- 
dered, albeit with very large intrinsic scales, the nematic 
order parameter Q shows a discontinuous jump at the 
transition. The disorder/polar transition occurs between 
the microscopically disordered phase D and the segre- 
gated phase P W avos ■ Like for the original Vicsek model, it 
is also discontinuous, as the polar order parameter takes 
finite, order 1 values as soon as the waves appear Q. 

We studied the two different P-N transitions present. 
The transition between the segregated phases P waV es and 
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FIG. 3: (color online) Direct transitions from nematic to polar 
order in the generalized Vicsek model, (a) transition between 
the segregated phases P wavcs and Nband at s — 0.1 occurring 
near r\ — 0.132 (time series of polar (top) and nematic (bot- 
tom) order parameter for L = 256). (b-d) transition between 
the homogenous phases Phom and Nhom at s = 0.01 occur- 
ring near r\ — 0.085. (b,c): order parameter curves at sizes 
L = 128, 256, 512, 1024. Panel (c) is a close-up of panel (b). 
(d): Binder cumulant curves at the same system sizes. The 
arrows indicate increasing system size. 



Nbands was studied as a function of 77 at s = 0.1. It 
is clearly discontinuous, as seen, e.g., in the characteris- 
tic flip-flop dynamics of the polar order parameter P in 
the transition region, leading to a bimodal distribution 
testifying of phase coexistence (Fig. [3^)- The transition 
between the homogeneous nematic and polar phases is 
more difficult to characterize. We performed a finite-size 
scaling study at s = 0.01, varying 77 around the transition 
point. The behavior of the P(rj) curves (Fig. \Sjp) indi- 
cates the premises of a discontinuous transition: although 
no discontinuity proper is present even at the largest size 
considered, these curves cross each other, suggesting that 
a jump may appear at still larger sizes. This conclusion 
is also borne out of the behavior of the Binder cumulant 
G = 1 - (P(i) 4 ) t /(3(P(i) 2 )?), which develops a deeper 
minimum as the system size is increased (Fig. [3ji) [13] ■ 

We now report on the properties of the metric-free, 
"topological" version of our generalized Vicsek model 
with competing ferromagnetic and nematic alignment. 
This case is of theoretical interest because no density- 
segregated phases are present in the purely ferromag- 
netic or nematic cases 17-19j], opening the way to con- 
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tinuous phase transitions (see below). It is also relevant 
in the context of collective motion. Common sense and 
some experimental/observational evidence indicate that 
in groups of higher organisms moving collectively (bird 
flocks, fish schools, crowds, etc.), one individual, rather 
than interacting with neighbors located within a given 
metric zone around itself, takes into account those in- 
dividuals forming some "angular landscape" to perform 
navigation decisions 
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It was suggested for instance 
that starlings interact with their 7-8 nearest neighbors, 
irrespective of the flock density . Careful study of fish 
trajectory data showed that in some schools the stimu- 
lus/response function of fish can be modeled by inter- 
actions with their first Voronoi neighbors (those whose 
polygons, in a Voronoi tessellation of space, form the 
first shell around the focal fish) ilfij. Here, following a 
study of a metric-free version of the Vicsek model 
we consider these Voronoi neighbors, and we use 1 



17], 



•vec- 



torial" noise, which means that Eq.([T]) is replaced by 

where £ is a ran- 
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arg 



domly oriented unit vector and TV} is the current number 
of (Voronoi) neighbors of particle j. 



As shown in [17|, |18| , global density drops out of the 
problem in metric- free models (it can be scaled out). In 
practice we worked at unit density, with vq = 0.5. The 
phase diagram of this metric-free version of our gener- 
alized Vicsek model is shown in Fig. 0^ for L = 128 
(diagrams obtained at larger sizes are nearly indistin- 
guishable). As expected, no density-segregated phases 
are present, leaving a diagram qualitatively similar to 
that of the equilibrium case. The P-D and N-D tran- 
sitions are now continuous, as expected from 
where the s = 1 and s = cases were studied. 
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A de- 
tailed study of the associated critical exponents will be 
presented elsewhere. (For the P-D transition, they are 
consistent with the values reported in [13] ■) We studied 
the direct P-N transition at s = 0.05 by finite-size scaling. 
The (polar) order parameter curves now do not cross each 
other (Fig. 0b), and the Binder cumulant curves show a 
minimum which, after deepening at small sizes, eventu- 
ally start receding at the largest sizes we could probe 
(Fig. Hfc). These qualitative facts point to a continuous 
transition. Quantitatively, our estimates of critical expo- 
nents are as follows: The crossings of the Binder cumu- 
lant curves (near G ~ | , not visible on Fig. |4}:) converge 
to an asymptotic threshold rj c — 0.4530(1) with an expo- 
nent 1/v = 1.00(5). The location of the maxima of the 
susceptibility also converge to the same estimated thresh- 
old with the same estimated in excellent agreement 
with the Ising value v — 1 (Fig. [4}i). The peak values 
of the susceptibility diverge with system size with expo- 
nent j/u = 1.74(2), the Ising value being | (not shown). 
The order parameter decreases algebraically at the esti- 
mated critical point with exponent (3/v = 0.126(3), in 
close agreement with the Ising value | (Fig. 0^) . Thus, 
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FIG. 4: (color online) Metric-free version of the model where 
particles interact with their Voronoi neighbors, (a): phase 
diagram obtained for L = 128. Legends as in Fig. [2kh. (b- 
d): finite-size scaling study of the P-N transition occuring at 
s = 0.05. (b) and (c): P(n) and G(rj) for L = 24, 32, 64, 128, 
and 256 (the arrows indicate increasing system size). The 
order parameter curves do not cross each other; the Binder 
cumulant curves show a minimum which first deepens then 
start receding as L increases. These curves do cross each 
other when G ~ | (not shown), (d): scaling of finite-size 
distance to asymptotic threshold Ari c (L) — \r)c{L) — n^°\ for 
estimated threshold = 0.4530(1). Red squares: r] c {L) is 
the location fo the maximum of susceptibility at size L. Black 
circles: r/ c (L) is the crossing point of two G(rf) curves at sizes 
Li and L2 with L = \/L\L2- The dashed blue line has slope 
-1. (e): scaling of order parameter at rj = r/^° . The dashed 
blue line has slope — |. 



in spite of rather strong finite-size effects as testified by 
the behavior of the Binder cumulant curves on the ne- 
matic side, all estimated critical exponents are very close 
to their Ising universality class values, as expected from 
studies of the equilibrium generalized XY model [9( . 

Before summarizing, we discuss briefly the derivation 
of a continuous theory describing active matter systems 
with competing ferromagnetic and nematic alignment in- 
teractions. Such a theory must a priori be in terms of a 
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polarity field P and a nematic tensorial field Q if it is to 
account for nematically-ordered phases. Baskaran and 
Marchetti have proposed rather complicated such equa- 
tions for the case of self-propelled rods interacting via 
steric exclusion . We have followed the somewhat sim- 
pler route of the "Boltzmann" approach used in 18l 21|. 
which is particularly adapted to dilute Vicsek-like mod- 
els, for the case of competing ferromagnetic and nematic 
alignment 22j ■ While details will be published elsewhere 
|23| . we only mention here that the obtained equations 
for P and Q fail to account for the phase diagrams pre- 
sented here. They confirm the existence of a nematic 
phase at finite s, but are not suitable to describe the po- 
lar phase, probably since the truncation scheme used is 
only valid near onset of nematic order. 

In conclusion, we have studied a Vicsek-style model 
with competing ferromagnetic and nematic alignment in 
both a metric and a "metric-free" version where interac- 
tions take place with Voronoi neighbors. We have shown 
that the fully nematic case of this out-of-equilibrium XY 
model ("self-propelled rods") resists some bias toward 
ferromagnetic alignment, thus conferring some robust- 
ness to its nematically-ordered phases and allaying our 
initial concern about friction effects inducing a "polar 
bias" in aligning collisions of elongated objects in ex- 
periments. In the metric case, the direct polar/nematic 
transition has been found discontinuous, in line with the 
order/disorder transitions, whereas the metric- free ver- 
sion exhibits an Ising-class continuous transition, as in 
equilibrium. We have signalled that a simple derivation 
of a continuous theory able to account for all observed 
facts is not easy, and constitutes an important future step 
for putting these results on firmer ground. 
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